Can f(R) Modified Gravity Theories Mimic a ACDM Cosmology? 
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We consider f(R) modified gravity theories in the metric variation formalism and attempt to 
reconstruct the function f(R) by demanding a background ACDM cosmology. In particular we 
impose the following requirements: a. A background cosmic history H(z) provided by the usual 
flat ACDM parametrization though the radiation (w e // = 1/3), matter (tu e // = 0) and deSitter 
(w e ff = —1) eras. b. Matter and radiation dominate during the 'matter' and 'radiation' eras 
respectively i.e. Q m = 1 when w e ff = and O r = 1 when w e ff = 1/3. We have found that 
the cosmological dynamical system constrained to obey the ACDM cosmic history has four critical 
points in each era which correspondingly lead to four forms of f(R). One of them is the usual 
general relativistic form f(R) = R — 2A. The other three forms in each era, reproduce the ACDM 
cosmic history but they do not satisfy requirement b. stated above. 

PACS numbers: 



I. INTRODUCTION 

There is accumulating observational evidence based 
mainly on Type la supernovae standard candles [l[ and 
also on standard rulers that the universe has entered 
a phase of accelerating expansion at a recent cosmolog- 
ical timescale. This expansion implies the existence of 
a repulsive factor on cosmological scales which counter- 
balances the attractive gravitational properties of matter 
on these scales. There have been several theoretical ap- 
proaches [J, [5[ towards the understanding of the origin of 
this factor. The simplest such approach assumes the ex- 
istence of a positive cosmological constant which is small 
enough to have started dominating the universe at re- 
cent times. The predicted cosmic expansion history in 
this case (assuming flatness) is 



H{zf 



where fio 



= Hi [fiom(l + zf + fior(l + zf + n A ] 

(1.1) 

£rs ^- ~ 10 4 is the present energy den- 



sity of radiation normalized over the critical density for 
flatness p cr it- Also fiom = ~ 0.3 is the normal- 

ized present matter density and Q\ = 1 — Qom — ^0r is 
the normalized energy density due to the cosmological 
constant. This model provides an excellent fit to the cos- 
mological observational data and has the additional 
bonus of simplicity and a single free parameter. Despite 
its simplicity and good fit to the data this model fails to 
explain why the cosmological constant is so unnaturally 
small as to come to dominate the universe at recent cos- 
mological times. This fine tuning problem is known as 
the coincidence problem. 

In an effort to address this problem two classes of mod- 
els have been proposed: The first class assumes that gen- 
eral relativity is a valid theory on cosmological scales 
and attributes the accelerating expansion to a dark en- 
ergy component which has repulsive gravitational prop- 
erties due to its negative pressure. The role of dark en- 



ergy is usually played by a minimally coupled to grav- 
ity scalar field called quintessence]^. Alternatively, the 
role of dark energy can be played by various perfect flu- 
ids (eg Chaplygin gas_0|), topological defects 0, holo- 
graphic dark energy etc. The second class of mod- 
els attributes the accelerating expansion to a modifica- 
tion of general relativity on cosmological scales which 
converts gravity to a repulsive interaction at late times 
and on cosmological scales. Examples of this class of 
models include scalar-tensor theories (Tol ITTI|, f(R) modi- 
fied gravity theories fl2j. braneworld models jl3| etc. An 
adva ntage o f models in this class is that they naturally 
allow [lOL for a superaccelerating expansion of the uni- 
verse where the effective dark energy equation of state 
w = £ crosses the phantom divide line w = — 1. Such a 
crossing is consistent with current cosmological data [HI] . 

Most of the models in both classes require the existence 
of arbitrary new degrees of freedom whose role is usually 
played by effective scalar fields. This is not a welcome 
feature because the degrees of freedom are to some extend 
arbitrary with respect to either their origin and/or their 
dynamical properties. Their predictive power is therefore 
usually dramatically diminished. 

A partial exception to this rule is provided by modified 
f(R) theories of gravity. In these theories the Ricci scalar 
R in the general relativistic Lagrangian is replaced by an 
arbitrary function f(R) leading to an action of the form 



S = 



-/(i?) + £ rad 



(1.2) 



where C m and £ ra d are the Lagrangian densities of matter 
and radiation and we have set 8irG = 1. These f(R) 
theories arise in a wide range of different frameworks: In 
quantum field theories in curved SDacetime[16||, in the low 
energy limit of the D = 10 superstring theory |17j . in the 
vacuum action for the Grand Unified Theories (GUTs) 
etc. 

It has been demonstrated [l8[ that for appropriate 
forms of f(R) the action (|1.2p can naturally produce 
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accelerating expansion at late times in accordance with 
Snla data[19|. The advantage of these theories is that no 
extra arbitrary degree of freedom is introduced and the 
accelerating expansion is produced by the Ricci scalar 
(dark gravity) whose physical origin is well understood. 
On the other hand, the main disadvantage of these the- 
ories is that (like most modified gravity theories) they 
are seriously constrained by local gra vity experiments 
HI HI HI- In fact it can be shown[2(l that f(R) mod- 
els are equivalent to scalar-tensor theories with vanish- 
ing Brans Dicke parameter (u> = 0) and a special type 
of potential. Since solar system tests of general relativ- 
ity imply u> > 4 x 10 4 [23|, these theories can only be 
consistent with observations if they are associated with 
a large (infinite) effective mass of the scalar R. It has 
been shown [24( that specific forms of the function f(R) 
can provide an infinite effective mass needed to satisfy 
solar system constraints and can also produce late time 
accelerating expansion. 

The reduction of f(R) theories to a special class of 
scalar-tensor theories implies that in principle the recon- 
struction of f(R) from a particular cosmic history H (z) 
can be performed in a similar way as in the case of scalar- 
tensor theories [Tol fl4|. However, the non-existence of a 
Brans Dicke parameter requires some modifications of the 
reconstruction methods especially when the reconstruc- 
tion extends through the whole cosmic history through 
the radiation and matter eras. The dynamical systems 
approach followed in the present study illustrates these 
modifications 

The construction of cosmologically viable models in- 
corporating late accelerating expansion based on f(R) 
theories has been an issue of interesting debate du ring 
the past year. This debate originated from Ref. [251 ] 
which demonstrated that f(R) theories that behave as a 
power of R at large or small R are not cosmologically vi- 
able because they have the wrong expansion rate during 
the matter dominated era (a ~ t 1 / 2 instead of a ~ i 2 / 3 ). 
This conclusion was challenged in Ref. [26[ which claimed 
that wide classes of f(R) gravity models including mat- 
ter and acceleration phases can be phenomenologically 
reconstructed by means of observational data. The de- 
bate continued with the recent Ref. [27| where a detailed 
and general dynamical analysis of the cosmological evo- 
lution of f(R) theories was performed. It was shown that 
even though most functional forms of f(R) are not cos- 
mologically viable due to the absence of the conventional 
matter era required by data, there are special forms of 
f(R) that can be viable (consistent with data) for appro- 
priate initial conditions. 

In the present study we perform a generic model in- 
dependent analysis of f(R) theories. Instead of speci- 
fying various forms of f(R) and finding the correspond- 
ing cosmological dynamics, we specify the cosmological 
dynamics to that of the ACDM cosmology and search 
for a possible corresponding form of f(R). We thus at- 
tempt to reconstruct f(R) from the background cosmo- 
logical dynamics. In particular we consider the general 



autonomous system for cosmological dynamics of f(R) 
theories and study the dynamics of f(R) using as input 
a ACDM cosmic expansion history. Our study is per- 
formed both analytically (using the critical points and 
their stability) and numerically by explicitly solving the 
dynamical system. The results of the two approaches 
are in good agreement since the numerical evolution of 
f(R) follows the evolution of the 'attractor' (stable crit- 
ical point) of the system for most initial conditions. As 
we point out in the next section however the physical 
significance of this 'attractor' should be interpreted with 
care since it is an artifact of the allowed perturbations in 
the form of the physical law f(R). 

The structure of this paper is the following: In the next 
section we derive the autonomous system for the cosmo- 
logical dynamics of f(R) theories. Using as input a par- 
ticular cosmic history H(z) (eg ACDM) we show how can 
this system be transformed so that its solution provides 
the dynamics and functional form of f(R). We also study 
the dynamics of this transformed system analytically by 
deriving its critical points and their stability during the 
three eras of the cosmic background history (radiation, 
matter and deSitter). We find that there are 'attractor' 
critical points for each era which allow an analytical pre- 
diction of the dynamics of the system. We also confirm 
this analytical prediction by a numerical solution of the 
system demonstrating that the evolution of the system 
is independent of the initial conditions. In section III we 
use the solution of the above system to reconstruct the 
cosmological evolution and functional form of the func- 
tion f(R). We also demonstrate the agreement between 
the analytical and numerical reconstruction of f(R). Fi- 
nally in section IV we conclude, summarize and refer to 
future prospects of this work. 

II. DYNAMICS OF f(R) COSMOLOGIES 

We consider the action (|1.2[) describing the dynamics 
of f(R) theories in the Jordan frame [12| . In the context 
of flat Friedman-Robertson- Walker (FRW) universes the 
metric is homogeneous and isotropic ie 

ds 2 = -dt 2 + a 2 (t)dx 2 (2.1) 

and variation of the action (|1.2p leads to the following 
dynamical equations which are the generalized Friedman 
equations 

5FH 2 = p m + p rad + l(FR-f)- 3HF (2.2) 

-2FH = p m + -Prad + F — HF (2.3) 

where F = ^ and p m , p ra d represent the matter and 
radiation energy densities which are conserved according 
to 

p m + 3H Pm = , (2.4) 
Prad + 4Hp rad = . (2.5) 
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In order to study the cosmological dynamics implied by 
equations (|2.2p , (|2.3[) we express them as an autonomous 
system [27I] of first order differential equations. To achieve 
this, we first write (12.211 in dimensionless form as 



I = 



Prad 



R 



I 



where 



3FH 2 3FH 2 6H 2 6FH 2 



1 d 



F' 

IF (2-6) 



dbxa dN H dt 
We now define the dimensionless variables x\, 



x\ 

X4 



EL 

' ~F 



6FH 2 



R H 1 
6H 2 ~ ~H + 



Prad 



3FH 2 

where in (|2. 10|) we have used the fact that 
R = 6 (2H 2 + h) = 6 (2H 2 + H'H) 



(2.7) 

(2.8) 
(2.9) 
(2.10) 
(2.11) 

(2-12) 



and we can associate X4 with Q r and x± + X2 + X3 = Qde 
with curvature dark energy (dark gravity) . Defining also 



3jrjp we can write equation (|2.6[) as 
O m = 1 - xi — x 2 — x 3 — X4, 



(2.13) 



We may now use (|2.7[) to express 



as 



IT 



1 4^ + ELL + ELEL 

2 \FH 2 3 FH 2 F H F 



-1 — X3 — 3x 2 + %1 + X4 



EL 

~F 



(2.14) 



(2.15) 



Also, differentiating X4 of (|2.1ip with respect to N we 
have 



Prad Prad F' 2p ra d H' 

3FH 2 _ 3FH 2 F ~ 3FH 2 H 



x'a = —2X3X4 + X1X4 



(2.16) 



(2.17) 



where we have made use of (|2.5[) . Similarly, differentiat- 
ing (|2.9p with respect to N we find 



c , = X1X3 _ X2 (2x 3 - Xl -4) (2.I8) 



where 



/' ~ f,R 



(2.19) 



and implies derivative with respect to R. Finally dif- 
ferentiating (|2.10p with respect to N we find 
, X1X3 



2x 3 (x 3 - 2) 



(2.20) 



The autonomous dynamical system (|2TT5j) . (|STT5]) . (l2~2T)f . 
(|2.17[) is the general dynamical system that describes the 
cosmological dynamics of f(R) theories. It has been ex- 
tensively studied in Ref. [27[ for various cases of f(R) 
(or equivalently various forms of m) and was found to 
lead to a dynamical evolution that in most cases is in- 
compatible with observations since it involves no proper 
matter era. Some forms of f(R) however were found to 
lead to a cosmological evolution that is potentially con- 
sistent with observations. In order to investigate such 
cases in more detail we follow a different approach. In- 
stead of investigating the above autonomous system for 
various different behaviors of m(/(i?)) we eliminate m 
from the system by assuming a particular form for H(N) 
(ie X3(N) (see f|2 . 10[) 'l ) consistent with cosmological ob- 
servations. Once #3 (AT) is known we can solve (|2.20p for 
^^p- and substituting in (|2.18j) we find 

x' 2 = -x' 3 - 2x 3 (x 3 - 2) - x 2 (2x 3 - Xl - 4) (2.21) 



which along with (|2.15|) and (|2.17p consist a new dynam- 
ical system which is independent of m. The study of this 
system will be our focus in what follows. 

The results of our analysis do not rely on the use of 
any particular form of 2:3 (N) (ie H(z)). They only re- 
quire that the universe goes through the radiation era 
(high redshifts), matter era (intermediate redshifts) and 
acceleration era (low redshifts). The corresponding total 
effective equation of state 



is 



2 H'{N) 



w e ff — — Radiation Era 

o 



(2.22) 



(2.23) 



w ef f = Matter Era 
w e ff = —1 deSitter Era 

For the sake of definiteness however, we will assume a 
specific form for H(z) corresponding to a ACDM cos- 
mology (|1.1|) which in terms of N takes the form 

H{N) 2 = Hi [fl 0m e- 3N + n 0r e- 4N + A ] (2.24) 

where N = lna = — ln(l + z) and Qa = 1 — ^om — Qo r - 
We can use (|2~T0|) and (|2~24"1) to find x 3 (N) as 



X3 (N)=2--- 



Clo m e 



-3N 



-47V 



2 n 0m e- 3N + n 0r e-* N + (1 - n Qm - n 0r ) 

(2.25) 

The crucial generic properties of X3 (N) are its values at 
the radiation, matter and deSitter eras: 



x 3 (N) =0 N < N, 
xs(N) = L 



N rm <N <N„ 



x 3 (N) = 2 N > N n 



(2.26) 
(2.27) 
(2.28) 
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where N rm ~ — ln^^ and N m \ ~ -|ln^- are the 
iV values for the radiation-matter and matter-deSitter 
transitions. For Slo m = 0.3, f2or = 10~ 4 we have N rm ~ 
—8, iV m A — —0.3. The transition between these eras 
is model dependent but rapid and it will not play an 
important role in our analysis. 

It is straightforward to study the dynamics of the 
system ([235]) . (j2~2"T1) . (|2~T7)) by setting x\ = to find 
the critical points and their stability in each one of the 
three eras corresponding to (|2.26|) - (|2.28[) . Notice that 
even though this dynamical system is not autonomous 
at all times it can be approximated as such during the 
radiation, matter and deSitter eras when X3 is approxi- 
mately constant. The critical points and their stability 
are shown in Table I. 

The stability analysis of Table I assumes that X3 — 
const and therefore it is not identical to the full stability 
analysis where X3 would be allowed to vary. The usual 
stability analysis of cosmological dynamical systems as- 
sumes a particular cosmological model (eg a form of f(R) 
or m) and in the context of this 'physical law', the sta- 
bility of cosmic histories H(N) is investigated. In this 
context clearly a stable cosmic history is the one pre- 
ferred by the model. 

In the reconstruction approach however the stability 
analysis has a very different meaning. Here we do not 
fix the model f(R) ('physical law'). Here we fix the cos- 
mic history and allow the physical law f(R) to vary in 
order to predict the required cosmic history. Thus our 
stability analysis concerns the 'physical law' f(R) and 
not the particular cosmic history. Since the physical law 
is usually fixed by Nature the instabilities we find are 
not physically relevant but they are only useful to under- 
stand analytically the phase space trajectories we obtain 
numerically. The physically interesting quantities are the 
values of the critical points we find in each era in the 
context of the ACDM cosmic history. These tell us the 
possible physical laws f(R) that can reproduce a ACDM 
cosmic history. As shown in Table I, one of these laws is 
clearly the general relativistic f(R) = R — 2A. 

The important point to observe in Table I is that in 
each era there are four critical points one of which corre- 
spond to the general relativistic f(R) = R — 2A. Three 
of the four critical points in each era are not stable. This 
however does not imply that these points are not cosmo- 
logically relevant. These instabilities are not instabilities 
of the trajectory H(N) (which we keep fixed) but of the 
forms of f(R) which is allowed to vary. Thus they are 
not so relevant physically since in a physical context f(R) 
is assumed to be fixed a priori. The 'attractor' critical 
points of Table I are relevant only for technical reasons 
since they allow a comparison between a numerical evo- 
lution and an analytical prediction of the evolution of the 
system. In a more realistic situation where the perturba- 
tions of the 'physical law' f(R) would be turned off, all 
critical points would correspond to valid reconstructions 
of a ACDAI cosmic history. 

If we allow for f(R) perturbations (but not of H(N) 



=? -0.4 



-0.8 



-30 -25 -20 -15 -10 -5 5 

N 

FIG. 1: The effective equation of state w e ff(N) imposed on 
the dynamical system (obtained from (12 . 22 fl using Q2.24j ). 
The long dashed red line starts deep in the radiation era while 
the short dashed blue line starts in the matter era and ignores 
radiation (fior = 0). 

perturbations), the evolution of the system is determined 
by just following the evolution of the 'attractors' of Table 
I through the three eras. This evolution is presented in 
Table II showing the 'attractors' in each era. We stress 
however that this is not necessarily a preferred cosmolog- 
ical trajectory for the reasons described above. 

The 'standard' critical points are shown in Table III 
and they are the only critical points that have in addition 
to the correct expansion rate properties, the required val- 
ues of (fi r , fi m , CIde) m each era. As shown in the next 
section these saddle critical points reconstruct the gen- 
eral relativistic f(R) ie f(R) = R — 2A. It is therefore 
clear that nonlinear f(R) theories can produce an ob- 
servationally acceptable cosmic history but not with the 
required values of (f2 r , fi m , SIde) in each era. We should 
stress that our analysis has not excluded the possibility 
of physical values of (f2o r , Slom, ^de) in the case of cos- 
mic histories oscillating around the anticipated w e ff in 
each era or a w e ff that is continuously evolving. These 
special cases however maybe severely constrained obser- 
vationally. 

To confirm the dynamical evolution implied by the 'at- 
tractors' of Table I, we have performed a numerical anal- 
ysis of the dynamical system (|2.15[) . (|2.2ip , (|2. 1T[) using 
the ansatz (|2.25[) for X3 with fiom = 0.3 and f^or = 10~ 4 . 
This ansatz for x$(N) leads to the w e /f(N) shown in Fig. 
1. We have set up the system initially, close to the 'stan- 
dard' radiation era saddle point (0, 0, 0, 1) and allowed it 
to evolve. As seen in Fig. 2 and Fig. 3 much before the 
onset of the matter era (N = N rr ~ — 25 < — 8 ~ N rm ) 
the slow (but non-zero) evolution of x$(N) forces the 
phase space trajectory to depart from the saddle point 
(0, 0, 0, 1) and head towards the radiation era stable 'at- 
tractor' (—4, 5, 0, 0) where it stays throughout the rest of 
the radiation era (w e ff — g). 
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TABLE I: The critical points of the system (2.15), (2.21), (2.17) and their stability in each one of the three eras. Stable 
points (attractors) have only negative eigenvalues, saddle points have mixed sign eigenvalues and unstable points have positive 
eigenvalues. 



Era 


N Range 


Xl 


X2 


X3 


x± 


Eigenvalues 






-1 











(3,-2,-1) 


R,3.di&t ion 


TV < -In » 0m 

"Or 


1 


o 


o 


o 


(5 2 1) 


W e ff = h 
j J 3 




-4 


5 








(-5,-4,-3) 















1 


(4,-1,1) 






1 
1 


-o/ o 


1 /9 




^4.ooD, 1,U. 1 14 ) 


Matter 







-1/2 


1/2 





(3.386,-1,-0.886) 


W ef f = 




0.886 


-0.386 


1/2 





(4.272,0.886,-0.114) 






-3.386 


3.886 


1/2 





(-4.386, -4.272, -3.386) 









-1 


2 





(-4, -3, 1) 


deSitter 




-1 





2 





(-5, -4, -1) 


Weff = -1 




3 





2 





(4, 3, -1) 






4 





2 


-5 


(5, 4, 1) 



TABLE II: The 'attractor' critical point in each era. 



Era 


N Range 




Xl 


X3 


X4 


Weff 






Qrad 


Radiation 


N < -ln^L 

"Or 


-4 


5 








1/3 





1 





Matter 


-ln^am- < N < -ln^- 

"Or "Om 


-3.386 


3.886 


1/2 











1 





deSitter 




-1 





2 





-1 





1 






TABLE III: The 'standard' saddle critical points in each era. These are also the points producing a linear general relativistic 
f(R) = R - 2A (see equation (EHS})- 



Era 


N Range 




Xl 


X3 


£4 


W eff 


"Om 


Qde 


^rad 


Radiation 


N < -ln^L 

"Or 











1 


1/3 








l 


Matter 


-ln^ < TV < -ln^ 





-1/2 


1/2 








1 








deSitter 


N > -iln^ 





-1 


2 





-1 





1 
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-1 

-2 
-3 
-4 



(a) 
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1 

-1 
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(b) 



/V 




FIG. 2: The evolution of the variables xi(N), X2(N), xs(N) and X4(N) for 'standard' radiation era initial conditions (dotted 
red line) and 'standard' matter era initial conditions (dashed blue line) . The perturbed trajectories are rapidly dragged by the 
stable 'attractors' of each era. The numerically obtained evolution along the 'standard' saddle points of Table III is also shown 
(continuous green line). The instabilities of this trajectory are bypassed by using the constrained system (|2.30[l - (|2.31|) instead 
of the full system Ip35) . (|2~2Tj) . (p^Ojl . 
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FIG. 3: The phase space trajectories on the Xi — X2 plane (Fig3a) and £3 — X4 plane (Fig3b) for 'standard' radiation era initial 
conditions (dotted red line) and 'standard' matter era initial conditions (dashed blue line). The trajectory corresponding to 
the numerically obtained evolution along the saddle points of Table III is also shown (continuous green line). The points Ai, 
A2, A3 correspond to the 'attractors' of each era (radiation, matter and deSitter respectively) while the points Bi, B2, B3 
correspond to the 'standard' critical points of each era (see Tables II and III). Notice that on the projection of Fig. 3b the 
'attractor' points A%, A3 coincide with the 'standard' critical points B2, B3. 
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Subsequently, when x^(N) enters the matter era 
{w e ff — 0) at N rm ~ —8, the trajectory follows the evo- 
lution of the 'attractor' fixed point and heads towards 
the matter era 'attractor' (—3.386,3.886,0.5,0) ignoring 
the saddle point (0, —1/2, 1/2,0) of the 'standard' matter 
era. Finally when the matter era is over, the trajectory 
heads towards the deSitter 'attractor' (—1,0,2,0) which 
is also distinct from the 'standard' deSitter saddle point 
(0, —1, 2, 0). Notice that the deSitter 'attractor' is incon- 
sistent with observations due to the implied large vari- 
ation of the effective Newton's constant G(N) = y^vj 
even though this inconsistency could be ameliorated by 
'chameleon' type mechanisms [2{|. The evolution of 
(f2 r , f2 m , Qde) corresponding to the phase space trajec- 
tories of Figs. 2 and 3 is shown in Fig. 4a. Notice that 
Q, m = throughout the 'attractor' evolution of the sys- 
tem and the w e ff = of the matter era is induced by 
curvature dark gravity excitations. 

We have also tested initial conditions in the matter era 
starting the evolution on the saddle point (0, — 5,5,0) 
corresponding to the 'standard' matter era. In this case 
we also ignore radiation setting Q 0r = 0. We get an evo- 
lution of the system (see Figs 2, 3, 4b) which stays on 
the 'standard' matter era fl m — 1 for about 3 expan- 
sion times but before the onset of the acceleration era it 
gets absorbed by the 'attractor' towards the nonstandard 
deSitter critical point (—1,0,2,0). 

The above evolution along the 'attractor' critical 
points is a result of the 'physical law' f{R) perturbations. 
We can also reproduce trajectories that go through criti- 
cal points that are not stable by turning off these pertur- 
bations. For example we can recover the saddle critical 
point sequence 



(0, 0, 0, 1) - (0, ~ -i 0) -> (0, -1, 2, 0) (2.29) 



by fixing x x = in the system (gHHJ) , (|2~2T|) . ([2~T7]) and 
reducing it to the system 



1 — X3 — 3x 2 + 







-2x^X4 = x'a 



(2.30) 
(2.31) 



which can be easily solved using the ansatz (|2.25[) 

(see Fig. 2 and Fig. 3). The alternative approach of 
solving the decoupled pair (|2.18|) , (|2.17[) does not lead 
to the correct result because the f(R) perturbations are 
not turned off and the constraint is not respected in this 
case. 

It is straightforward to reconstruct the functions f(R) 
that correspond to the saddle general relativistic trajec- 
tory (|2.29p and to the 'attractor' sequence of Table II. 
The functional forms of f(R) may also be reconstructed 
on any one of the critical points of Table I. These tasks 
are undertaken in the next section. 



III. RECONSTRUCTION OF f(R) 

We now reconstruct the form of the function f(R) that 
corresponds to each one of the critical points of the sys- 
tem shown in Table I. This reconstruction is effectively an 
approximation of f(R) in the neighborhood of each crit- 
ical point. It is particularly useful because most of the 
dynamical evolution takes place close to the fixed points. 
Consider a critical point of the form (x%, £2, %3t X4). Us- 
ing (|2T8|) we find 



F = Fne- SlN 



(3.1) 



where Fo is a constant. We may eliminate N in favor 
of R using the input form of H(N) (equation (|2.24[) in 
(|2"?T2"|l to obtain (setting H§ = 1) 



R(N) = 3 [4ft A + f> 0m e~ 3Ar ] 



which leads to 



F = F 



iZ-12fi/ 

3f2om 



(3.2) 



(3.3) 



and by integration we get 

ffm 3F (3r> 0m )-»(i?-12Q A )»+ 1 

J\R) = - q r Jo (3.4) 

X\ -+- o 

where fo is an integration constant. Expressing (|3.4p in 
terms of N using (|3.2p we obtain 



f(N) = ^ oU °f, n +/ 



X\ 



(3.5) 



It is now straightforward to use the expressions for f(N) 
R(N) and H(N) to find x 2 (N) (equation (gU)), x 3 (N) 
(equation (|Q5]) j and x 4 (N) (equation (|2~TT|) ). We thus 
find 



and 



x 2 (N) 



x 4 {N) 



3n 0m 

2(2i+3)' 



-3JV 



+ 6F e J 



Ome 



-3N 



+ n 0r e- 4N + n A 



(Si-4)JV 



,-3N 



-AN 



(3.6) 



(3.7) 



while X3(N) is given by (|2.25|) . Using equations (|3.6p . 
(|2.25p and (|3.7p we may verify the x%, S3, X4 values of 
each critical point by considering the appropriate range 
of N in each era and the corresponding value of x±. By 
demanding consistency with the values of Table I we may 
obtain the values of the constants fo and Fo. 

As an example let's consider the sequence (|2.29p corre- 
sponding to the 'standard' cosmological eras (Table III). 
It is easy to see, using x± — and the appropriate range 
of N in ((331) . and Ij37f| that we obtain the correct 

values for x 2 , X3, X4 in the radiation and matter eras 
for any value of Fq, fo. In the deSitter era (N >> 1) 
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FIG. 4: The evolution of Q m (long-dashed), Q ra d (dotted) and Qde (short-dashed) parameters for 'standard' radiation era 
initial conditions (Fig4a) and 'standard' matter era initial conditions (Fig4b). In both cases the 'attractor' critical points of 
Table II rapidly take over and drag the trajectories. 



the value of fo needs to be fixed to get agreement with 
X2 = —1 of Table I. In particular from (|3.6j) we find 



X2 



fo 



6F n A 



= -1 



which implies 



fo = 6F Cl A 



(3.8) 



(3.9) 



Using now (|3.9[) and setting x± = in (|3.4p we recon- 
struct the expected result 



f(R)=F {R-6Q / 



(3.10) 



which is valid for all three eras since in this sequence the 
value of x\ remains constant. In a similar way we may 
reconstruct f{R) for any critical point in one of the three 
eras. 

We have therefore extended previous studies showing 
that f(R) theories can only be viable in very restricted 
cases by showing that even these restricted cases can not 
reproduce a viable ACDM cosmology where w e f / is con- 
stant during the matter and radiation eras and tto r , &m 
take their cosmologically anticipated values. It therefore 
becomes clear that if the accelerating expansion of the 
universe is due to physics in the gravitational sector it 
may probably have to be a more general theory than f(R) 
modified gravity. Such a theory could very well be scalar- 
tensor gravity (or equivalently coupled dark energy [30l]) 
whose cosmological dynamical properties and constraints 
need to investigated in detail. 

In the case of sequence of transitions among critical 
points which involve different values of x\ the reconstruc- 
tion can be done by either numerical determination of 
x\ (N) or by approximating it as a sequence of step func- 
tions. For example, the steps involved in the reconstruc- 
tion of the 'attractor' trajectory shown in Fig. 2 and Fig. 
3 are the following: 



1. Use (|2.8[) along with the numerical solution x\(N) 
to find the function F(N) = f, R (N) as 

F(N) = F e-^« Xl{N ' )dN ' (3.11) 

The numerical solution xi(N) of Fig. 2a can be 
approximated as a piecewise constant function with 
values determined by the corresponding 'attractors' 
of each cosmological era and by the initial condi- 
tions ie 



xi(N) 
xi(N) 
xi(N) 
xi(N) 



- 30 < N < N„ 
-4 N rr < N < N r 



(3.12) 
(3.13) 



-3.386 N rm < N < N mA (3.14) 
-1 N mA < N (3.15) 



(where N rr ~ —25) thus leading to an analytical 
approximation for F(N). The resulting form of 
\n(F(N)) in both the numerical reconstruction and 
its analytical approximation is shown in Fig. 5 (5a: 
'standard' radiation era initial condition, 5b: 'stan- 
dard' matter era initial condition). 



2. Use equation 
f(N) 



to find f(N) from F(N) ie 
6x 2 (N)F(N)H(N) 2 (3.16) 



where H(N) is given by (|2.24p . X2(N) is numeri- 
cally obtained and shown in Fig. 2b and F(N) is 
obtained in the previous step. As in the case of 
Xi(N), X2^N) can be analytically approximated as 

x 2 {N) =0 -30<A^<^ rr (3.17) 

x 2 {N) =5 N rr < N < N rm (3.18) 

x 2 (N) = 3.886 N rm < N < N mA (3.19) 

x 2 (N) ~ N mA < N (3.20) 

using the corresponding 'attractors' to obtain an 
analytical expression for f(N). 
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FIG. 5: The form of log(F(N)) in the numerical reconstruction (dashed lines) and its analytical approximation using both the 
interpolating steps 1-3 (dotted lines) and the application of the analytical expression (|3.5[1 valid in each era (thick green lines). 
The agreement between the three approaches is very good. 5a: 'standard' radiation era initial condition, 5b: 'standard' matter 
era initial condition. 
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FIG. 6: The form of f(R) in the numerical reconstruction (dashed lines) and its analytical approximation using both the 
interpolating steps 1-3 (dotted lines) and the application of the analytical expression (|3.4|l valid in each era (continuous green 
lines). 6a: 'standard' radiation era initial condition with the continuous green lines corresponding to radiation era (larger R) 
and matter era (smaller R). The deSitter era is not shown since it corresponds to a single point (constant R). 6b: 'standard' 
matter era initial condition with the continuous green line corresponding to matter era. The deSitter era is not shown since it 
corresponds to a point (constant R). 



3. The resulting form of f(N) can then be combined 
with equation (|3.2p for R(N) to reconstruct the 
function f(R). The resulting form of f(R) is shown 
in Fig. 6 for both the numerical reconstruction and 
its analytical approximation (6a: radiation era ini- 
tial conditions, 6b: matter era initial conditions). 

We can fit the reconstructed f(R) of Figs. 6a and 6b to 
the analytic form of equation (|3.4[) for each era respec- 
tively so as to find the parameters Fq, x± and /o- The 
results are shown in Table IV. 

Notice that the best fit values of X\ coincide with the 
corresponding 'attractor' critical points of Table II as ex- 
pected. This verifies the validity of the reconstructed 
f(R) expression from (|3.4| . A similar reconstruction 



TABLE IV: The parameters Fq, x\ and fo for the recon- 
structed f(R) of Figs. 6a and 6b . 



Radiation Era cond. 




Log(F ) Log(-f ) 


Radiation Era 


-3.99 


96.73 


90.92 


Matter Era 


-3.45 


92.17 


90.59 


deSitter era 


-1.01 


92.62 


94.03 


Matter Era cond. 


Xl 


Fo 


fo 


Matter Era 





1 


-1.51 


deSitter era 


-1.02 


4.83 


-16.56 
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analysis can be performed for any other sequence of criti- 
cal points. As discussed in section II any such sequence is 
equally interesting cosmologically since the existence of 
the 'attractor' is an artifact of the f(R) perturbations. 



IV. CONCLUSION-OUTLOOK 

We have shown analytically and numerically that non- 
linear /(i?) gravity theories can reproduce the back- 
ground expansion history H{z) indicated by observations 
even when f(R) does not reduce to general relativity at 
early times. In that case the universe gets dominated 
by dark gravity during its evolution as opposed to radia- 
tion or matter. This result relies on the values of all the 
critical points we found assuming only that the radiation 
era corresponds to a constant effective equation of state 
parameter w e ff — -| while for the matter era we have 

W eff = 0. 

Our analysis indicates f(R) models can be viable if 
f(R) deviates from general relativity at early times. 
Thus a viable f(R) theory must satisfy one of the fol- 



lowing: 

• Either f(R) reduces to general relativity at early 
times, but departs from general relativity at late 
times (a well known case|12j). 

• Or dark gravity in the forms derived in our paper 
mimics radiation or matter at both the background 
level and the perturbative level. The later would 
clearly require a separate analysis of perturbations 
of the model. 

Numerical Analysis: The mathematica files with 
the numerical analysis of this study may be found 
at http://leandros.physics.uoi.gr/frlcdm/frlcdm.htm or 
may be sent by e-mail upon request. 
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